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DAMA is searching for rare processes by developing and using several kinds 
of radiopure scintillators: in particular, Nal(Tl), liquid Xenon and CaF 2 (Eu). 
The main results are here summarized with particular attention to the inves¬ 
tigation of the WIMP annual modulation signature. 
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1 Introduction 


DAMA is devoted to the search for rare processes by developing and using low radioac¬ 
tivity scintillators. Its main aim is the search for relic particles (WIMPs: Weakly In¬ 
teracting Massive Particles). In addition, due to the radiopurity of the used detec¬ 
tors and of the installations, several searches for other possible rare processes are also 
carried out, such as e.g. on exotics, on ( 3(3 processes, on charge-non-conserving pro¬ 
cesses, on Pauli exclusion principle violating processes, on nucleon instability and on 
solar axionsg, g, g, g, g. g. g. g. 0, [H HI EDI ■ 

The main experimental set-ups running at present are: the ~ 100 kg Nal(Tl) set-up, 
the ~ 6.5 kg liquid Xenon (LXe) set-up and the so-called “R&D” apparatus. Moreover, 
a low-background germanium detector is operative underground for measurements and 
selections of sample materials. 

In the following the most recent results achieved with the ~ 6.5 kg LXe set-up and 
with the “R&D” apparatus will be briefly recalled, while a more dedicated discussion will 
be devoted to the present results on the investigation of the WIMP annual modulation 
signature by the ~ 100 kg Nal(Tl) set-up. 


2 The LXe experiment 


We pointed out the interest in using liquid Xenon as target-detector for particle dark 
matter search deep underground since ref. [|Hj . Several prototypes were built and related 
results published [lg. The final choice was to realize a pure liquid Xe scintillator directly 
collecting the emitted UV light and filled with Kr-free Xenon gas isotopically enriched. 
The detailed description of this set-up and of its performances have been described in 
ref. |U|. Kr-free Xenon enriched in 129 Xe at 99.5% has been used since time, while more 
recently the set-up has been modified in order to run alternatively either with this gas or 
with Kr-free Xenon enriched in 136 Xe at 68.8%. 

After preliminary measurements both on elastic and inelastic WIMP- 129 Xe scattering 
17], [18|, the recoil/electron light ratio and the pulse shape discrimination capability in 


a similar pure LXe scintillator have been measured both with Arn-B neutron source and 
with 14 MeV neutron generator [|lg. After some upgrading of the set-up, new results on 

the pulse shape 
|1]) . Moreover, in 


the WIMP search have been obtained gU|. In particular, in ref. 
discrimination in pure LXe scintillators has been exploited (see Fig. 

2000/2001 further measurements on the recoil/electron light ratio with 2.5 MeV neutron 
generator have been carried out at ENEA-Frascati; see ref. [gl] for details and compar¬ 
isons. Fig. U summarizes the measured values. The inelastic excitation of 129 Xe by Dark 
Matter particles with spin-dependent coupling has also been searched for preliminarily in 


ref. [|18] and, more recently, in ref. [g(J . 


Other rare processes have been investigated by filling the detector with the Kr-free 
Xenon gas enriched in 129 Xe at 99.5%. 
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Figure 1: Low energy distribution (statistics of 1763.2 kg x day ) measured with the vessel 
filled with 129 Xe; the continuous line represents the upper limits at 90% C.L. obtained for 


the recoil fractions 


191 


In particular, as regards the electron stability, limits on the lifetime of the electron 
decay in both the disappearance and the v e + 7 channels were set in ref. [|TJ. The limits 
on the lifetime of the latter channel has been more recently improved to: 2.0(3.4) • 10 26 
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Figure 2: Measured behaviour of the recoil/electron light ratio with recoil energy. Note 
that the energy of the data point of A 111 B neutron source is an average which refers also 
to much lower energies than those explored with the 2.5 MeV neutrons. 


2 












y at 90% (68%) C.L. [[[(J. Furthermore, new lifetime limits on the charge non-conserving 


electron capture with excitation of 129 Xe nuclear levels have been established to be in the 
range (1 — 4) ■ 10 24 y at 90% C.L. for the different excited levels of 129 Xe ||. The most 
stringent restrictions on the relative strengths of charge non-conserving processes have 
been derived: e\ r < 2.2 x 10 -26 and e 2 < 1.3 x 10~ 42 at 90% C.L. 


W - ---- — - 7 

Moreover, the nucleon and di-nucleon decay into invisible channels 


111 has been in¬ 


vestigated by searching for the radioactive daughter nuclei, created after the nucleon or 
di-nucleon disappearance in the parent nuclei. This new approach has the advantage of 
a branching ratio close to 1 and - if the parent and daughter nuclei are located in the 
detector itself - also an efficiency close to 1. The obtained limits at 90% C.L. are: r(p —> 
invisible channel) > 1.9 • 10 24 y; r(pp —» invisible channel) > 5.5 • 10 23 y and r(nn —> 
invisible channel) > 1.2 • 10 25 y. These limits are similar or better than those previously 
available; the limits for the di-nucleon decay in u T u T have been set for the first time; 


moreover, the limits are valid for every possible disappearance channel [11 


Finally as mentioned above, more recently the set-up has been modified to allow the 
use of Kr-free Xenon enriched in 136 Xe at 68.8%. Preliminary measurements have been 
carried out during 6843.8 hours |T3[ and new experimental limits have been obtained 


for (3(3 decay processes in 136 Xe, improving those previously available by factors ranging 
between 1.5 and 65. In particular, for the 3 possible channels without neutrinos the 
following half life limits (90% C.L.) have been achieved: 7.0 • 10 23 y for the channel (3(3 
0z/(0 + —> 0 + ); 4.2 • 10 23 y for the channel (3/3 0z/(0 + —■> 2 + ); 8.9 • 10 22 y for the channel 
(3(3 0z/M(0 + —> 0 + ). For comparison, we note that the obtained experimental limits on 
the half life of the process (3(3 0z/(0 + —> 0 + ) is lower only than the one obtained for the 
case of 76 Ge. Moreover, the limits (90% C.L.) on the channels (3(3 2z/(0 + —» 0 + ) 1.1 • 10 22 
y and (3(3 0z/M(0 + —> 0 + ) are at present the most stringent ones not only for the 136 Xe 
isotopes, but also for every kind of nucleus investigated so far either by active or by passive 
source method. Furthermore, upper bounds on the effective neutrino mass have been set 
considering various theoretical models for the evaluation of the elements of the nuclear 
matrix; they vary between 1.5 eV and 2.2 eV (90% C.L.). Finally, in the framework of 
the same models we have also obtained upper limits on the effective coupling constant 
Majoron - neutrino; they range in the interval: 4.8 -10~ 5 and 7.1 -10 -5 (90% C.L.). 

In conclusion, competitive results have been achieved by the DAMA LXe set-up by 
using Kr-free isotopically enriched Xenon gas. Further upgradings to improve the detec¬ 
tor performance are under consideration, while the data taking is continuing with both 
enrichments. 


3 Results from the “R&D” set-up 

The “R&D” set-up is used to measure the performances of prototype detectors and to 
perform small scale experiments. In particular, several experiments on the investigation 
of (3(3 decay processes in 136 Ce, 142 Ce, 40 Ca, 46 Ca, 48 Ca, 130 Ba, 106 Cd 0, |], 0 0 ( see Fig. 
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[D; other measurements are in progress and/or in preparation. 
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Figure 3: Light gray: experimental limits on the half-life of several /3/3-decay-candidate 
isotopes investigated by DAMA (90% C.L. with exception of the results on Barium and 
Cerium isotopes, which are at 68% C.L.) EMM- Dark gray: previous best experimental 
limits (90% C.L.). 


4 Results from the ~ 100 kg Nal(Tl) set-up 


The ~ 100 kg Nal(Tl) DAMA set-up allows us to explore several kinds of rare processes 
regarding WIMPs, processes forbidden by Pauli exclusion principle, exotics, charge non 
conserving processes and solar axions|| 0. H, |l^, [22|. |23], ^4| |25|, [26|. HU- HU |29H . A full 


description of the apparatus and its performances can be found in ref. . 

Let us preliminarily briefly summarize recent results - not referred to the WIMP search 
- obtained by means of this set-up. 

A search for processes normally forbidden by the Pauli exclusion principle was carried 
out and an improved limit for spontaneous emission rate of protons in 23 Na and 12 ‘I has 
been determined. In particular, the upper limit on the unit time probability of non-paulian 
emission of protons with energy E p > 10 MeV is 4.6 x 10 -33 s _1 and the corresponding 
lower limit on the mean life results 0.7 x 10 25 y for 23 Na and 0.9 x 10 25 y for 12, I [3]. 
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The same set-up has allowed to set a new limit on the electron stability by considering 
for the first time the ’’disappearance” channel of L-shcll electrons in Iodine. With a 
statistics of 19511 kg x day we have obtained: r > 4.2 x 10 24 y (68% C.L.) for the process 
e~ —> 3 u, major cm + v or anything invisible ||. 

Furthermore, new lifetime limits have been obtained on the charge non-conserving 
electron capture with excitation of 127 I and 23 Na nuclear levels. By analysing a statistics 
of 34866 kg x day, limits in the range: (1.5 — 2.4) x 10 23 y (up to two orders of magnitude 
higher than the ones previously available) have been obtained; also stringent restrictions 
on the relative strengths of the charge non-conserving (CNC) processes have been set [7j. 

Competing results have also been achieved by searching for neutral SIMPs: neutral 
particles with masses between few GeV and the GUT scale and cross sections on protons 
up to 10 -22 cm 2 , embedded in the galactic halo (/3 ~ 1CU 3 ). Although their large cross 
section, the neutral SIMPs interaction rate should be relatively low because of their low 
density in the galactic halo. Two ’’planes” of Nal(Tl) in the ~ 100 kg Nal(Tl) DAMA 
set-up have been used for this search, excluding neutral SIMPs with masses up to 4 x 10 16 
GeV |24|]. With the same set-up also a new limit has been set on the neutral nuclearites 
(a possible new form of matter containing roughly equal number of up, down and strange 
quarks.). The model independent upper limit at 90% C.L. on the neutral nuclearites flux 
is: $ < 1.9 x 10 -11 s -1 cm -2 sr -1 |24j. 

Recently this set-up has also allowed to set a new limit on solar axions by exploiting 
the Primakoff effect in Nal crystals. Analysing a statistics of 53437 kg x day the limit 
g ail < 1.7 x 10~ 9 GeV^ 1 (90% C.L.) [T2j has been obtained (see Fig. |]); the region 
with axion mass greater than 0.03 eV (not accessible by other direct methods) has been 
explored and KSVZ axions with mass greater than 4.6 eV have been excluded at 90 % 
C.L. (see Fig. [|). 

As regards the WIMP search, the recoil/electron light ratio and the pulse shape dis¬ 
crimination capability have been measured by neutron source and competitive results 
have been achieved on the WIMP-nucleus elastic scattering by exploiting the pulse shape 
discrimination technique in ref. |22|. Moreover, studies on the possible diurnal variation 


of the low energy rate in the data collected by the ~ 100 kg Nal(Tl) have been carried out. 
We recall that a diurnal variation of the low energy rate in WIMP direct searches can be 
expected because of the Earth’s daily rotation. In fact, during the sidereal day the Earth 
shields a given detector with a variable thickness, eclipsing the “wind” of Dark Matter 
particles but only in case of high cross section candidates (to which small halo fraction 
would correspond). By analyzing a statistics of 14962 kg x day no evidence of diurnal rate 
variation with sidereal time has been observed |23|; this result supports that the effect 


pointed out by the studies on the WIMP annual modulation signature (see later) would 
account for a halo fraction > 10 3 . 

Let us note that the highly radiopure ~ 100 kg Nal(Tl) set-up has been originally 
developed to investigate the WIMP annual modulation signature, thus the remaining 
part of this paper will be devoted to summarize the present results on this topics. 
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Figure 4: Exclusion 
The limit achieved 
and previous direct 


plot for the axion coupling constant, g a77 , versus the axion mass, m a . 
by this experiment is shown together with theoretical expectations 


searches for solar axions 12 


5 Investigation of the WIMP annual modulation sig¬ 
nature 

The WIMPs are embedded in the galactic halo; thus our solar system, which is moving 
with respect to the galactic system, is continuously hit by a WIMP “wind”. The quanti¬ 
tative study of such “wind” allows both to obtain information on the Universe evolution 
and to investigate Physics beyond the Standard Model. The WIMPs are mainly searched 
for by elastic scattering on target nuclei, which constitute a scintillation detector. In 


particular, the ~ 100 kg Nal(Tl) set-up |25|] has been realized to investigate the so-called 
WIMP “annual modulation signature”. In fact, since the Earth rotates around the Sun, 
which is moving with respect to the galactic system, it would be crossed by a larger WIMP 
flux in June (when its rotational velocity is summed to the one of the solar system with 
respect to the Galaxy) and by a smaller one in December (when the two velocities are 
subtracted). The fractional difference between the maximum and the minimum of the 
rate is expected to be of order of ~ 7%. The ~ 100 kg highly radiopure Nal(Tl) DAMA 
set-up [^] can effectively exploit such a signature because of its well known technology, of 
its high intrinsic radiopurity, of its mass, of the deep underground experimental site and 
of its suitable control of the operational parameters. The annual modulation signature is 
very distinctive as we have already pointed out [[Z|, ^7], ^3], [J(], [TT], [T^j . In fact, a 

WIMP-induced seasonal effect must simultaneously satisfy all the following requirements: 
the rate must contain a component modulated according to a cosine function (1) with one 
year period (2) and a phase that peaks around ~ 2 nd June (3); this modulation must be 
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found in a well-defined low energy range, where WIMP induced recoils can be present (4); 
it must apply to those events in which just one detector of many actually ”fires”, since 
the WIMP multi-scattering probability is negligible (5); the modulation amplitude in the 
region of maximal sensitivity must be <7% (6). Only systematic effects able to fulfil these 
6 requirements could fake this signature; therefore for some other effect to mimic such a 
signal is highly unlikely. 

Here the results obtained by analysing the data of four annual cycles investigating the 
annual modulation signature by means of the ~ 100 kg Nal(Tl) set-up are summarized 
(see Table []]). Cumulative analyses of a total statistics of 57986 kg x day have been carried 


Table 1: Released data sets 22, 26 


29, 31, 32 


period 

statistics {kg x day) 

DAMA/NAI-1 

4549 

DAMA/NaI-2 

14962 

DAMA/NaI-3 

22455 

DAMA/NaI-4 

16020 

Total statistics 

57986 

+ DAMA/NaI-0 upper limit on recoils by PSD 


out [^9|, |3(J, |3lJ properly including the physical constraint which arises from the measured 
upper limit on the recoil rate |22], |3(J. Both model independent and model dependent 
analyses have been performed. We remark that various uncertainties exist in every model 
dependent calculations concerning WIMP direct detection search (the same is for exclusion 
plot as well as for allowed regions in the WIMP-proton cross section versus WIMP mass 
plane). They are present on some general features such as e.g. the real behaviour of the 
WIMP velocity distribution (which we have considered so far isothermal and Maxwellian), 
on the adopted particle model and on the values taken for all the parameters needed in 
the calculations. 


6 Annual modulation signature and related studies 

6.1 Model independent analysis 

A model independent analysis of the data of the four annual cycles offers an immediate 
evidence of the presence of an annual modulation of the rate of the single hit events in 
the lowest energy interval (2-6 keV) as shown in Fig. |5|. There each data point has been 
obtained from the raw rate measured in the corresponding time interval, after subtracting 
the constant part. 
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Figure 5: Model independent residual rate for single hit events, in the 2-6 keV cumulative 
energy interval, as a function of the time elapsed since January 1-st of the first year of 
data taking. The expected behaviour of a WIMP signal is a cosine function with minimum 
roughly at the dashed vertical lines and with maximum roughly at the dotted ones. 



The x 2 test on the data of Fig. j| disfavors the hypothesis of unmodulated behaviour 
(probability: 4 • 10 -4 ), while fitting these residuals with the function A- cosu;(f — to), one 
gets: i) for the period T = 77 = (1-00 ± 0.01) year when t 0 is fixed at the 152.5 th day of 
the year (corresponding to ~ 2 June); ii) for the phase to = (144 ± 13) days, when T is 
fixed at 1 year. In the two cases A is: (0.022 ± 0.005) cpd/kg/keV and (0.023 ± 0.005) 
cpd/kg/keV, respectively. 

We have extensively discussed the results of the investigations of the known sources of 
possible systematics when releasing the data of each annual cycle; moreover, a dedicated 
has been released on possible systematics, where in particular the data of the 


paper 


30 


DAMA/NaI-3 and DAMA/NaI-4 running periods have been considered in quantitative 
evaluations. No known systematic effect or side reaction able to mimic a WIMP induced 
effect has been found as discussed in details in ref. [|30 . 


In conclusion, a WIMP contribution to the measured rate is candidate by the result of 
the model independent approach independently on the nature and coupling with ordinary 
matter of the possible WIMP particle. 


6.2 Model dependent analyses 

To investigate the nature and coupling with ordinary matter of a possible candidate, a 
suitable energy and time correlation analysis is necessary as well as a complete model 
framework. We remark that a model framework is identified not only by general astro- 
























physical, nuclear and particle physics assumptions, but also by the set of values used for 
all the parameters needed in the model itself and in related quantities (for example WIMP 
local velocity, vq, form factor parameters, etc.). 

At present the lightest supersymmetric particle named neutralino is considered the best 
candidate for WIMP. In supersymmetric theories both the squark and the Higgs bosons 
exchanges give contribution to the coherent (SI) part of the neutralino cross section, while 
the squark and the Z° exchanges give contribution to the spin dependent (SD) one. Note, 
in particular, that the results of the data analyses [^, |3l], [32| summarized here and in 
the following hold for the neutralino, but are not restricted only to this candidate. 

The differential energy distribution of the recoil nuclei in WIMP-nucleus elastic scat¬ 
tering can be calculated |22], |33| by means of the differential cross section of the WIMP- 
nucleus elastic processes 


da 


dE 


R 


( V , Er) = 


da 


dE 


+ 


da 


R . 


SI 


dE 


R, 


( 1 ) 


SD 


where v is the WIMP velocity in the laboratory frame and Er is the recoil energy. 

In the following, the results obtained by analysing the data in some of the possible 
model frameworks are summarized. 

I. WIMPs with dominant SI interaction in a given model framework 

Often the spin-independent interaction with ordinary matter is assumed to be domi¬ 
nant since e.g. most of the used target-nuclei are practically not sensitive to SD interac¬ 
tions as on the contrary 23 Na and 12 'I are. Therefore, first model dependent analyses of the 
data have been performed by considering a candidate in this scenario, that is neglecting 

the term (m) SD in eq- (!)■ 

A full energy and time correlation analysis - properly accounting for the physical con¬ 
straint arising from the measured upper limit on recoils [p3, |30| - has been carried out 


in the framework of a given model for spin-independent coupled candidates with mass 
above 30 GeV. A standard maximum likelihood method has been used]]. Following the 
usual procedure we have built the y log-likelihood function, which depends on the exper¬ 
imental data and on the theoretical expectations; then, y is minimized and parameters’ 
regions allowed at given confidence level are derived. Note that different model frame¬ 
works (see above) vary the expectations and, therefore, the cross section and mass values 
corresponding to the y minimum, that is also the allowed region at given C.L.. In partic¬ 
ular, the inclusion of the uncertainties associated to the models and to every parameter 
in the models themselves as well as other possible scenarios enlarges the allowed region as 
discussed e.g. in ref. 123] for the particular case of the astrophysical velocities. Also in the 


case considered here the minimization procedure has been repeated by varying the WIMP 
local velocity, Vq, from 170 krn/s to 270 krn/s to account for its present uncertainty. For 


* Substantially the same results are obtained with other analysis approaches such as e.g. the Feldman 
and Cousins one. 
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example, in the model framework considered in ref. [^9], the values rriw = (72l}|) GeV 
and £(Jsi — (5-7 ± 1.1) • 10 -6 pb correspond to the position of y minimum when Vq = 
170 krn/s, while rriw = (43±g 2 ) GeV and C,asi = (5.4 ± 1.0) ■ 10~ 6 pb when vq = 220 
krn/s. Here, £ is the WIMP local density in 0.3 GeV cm -3 unit, &si is the point-like 
SI WIMP-nucleon generalized cross section and rriw is the WIMP mass. Fig. |6] shows 
the regions allowed at 3cr C.L. in such a model framework, when the uncertainty on Vq 
is taken into account (solid contour) and when possible bulk halo rotation is considered 
(dashed contour). For simplicity, no other uncertainties on the used parameters have been 
considered there (some of them have been included in the approach summarized in the 
next subsection [[Jl|]), which obviously further enlarge the allowed region with the respect 
both to £a si and rriw 


_Q 

Q_ 


b 

SUM 


10 " 


10 " 


10 


-6 


10 


-7 



10 


20 


50 


100 

(GeV) 


200 


500 


1000 


Figure 6: Regions allowed at 3<r C.L. on the plane £asi versus m. w for a WIMP with 
dominant SI interaction and mass above 30 GeV in the model framework considered in 
ref. |29|: i) when u 0 uncertainty (170 km/s < Vq < 270 km/s; continuous contour) 
has been included; ii) when also a possible bulk halo rotation as in ref. (dashed 

contour) is considered. As widely known, the inclusion of present uncertainties on some 
other astrophysical, nuclear and particle physics parameters would enlarge these regions 
(varying consequently the position of the minimum for the y log-likelihood function). 


A quantitative comparison between the results of the model independent analysis and 
of this model dependent analysis has been discussed in ref. . 

In conclusion, the observed effect investigated in terms of a WIMP candidate with 
dominant SI interaction and mass above 30 GeV in the model framework considered in 
ref. [|9|, supports allowed WIMP masses up to 130 GeV (1 a C.L.) and even up to 
180 GeV (1 a C.L.) if possible dark halo rotation is included, while lower £asi would be 
implied by the inclusion of known uncertainties on parameters (for example, on the form 
factor parameters) and model features. 

Theoretical implications of these results in terms of a neutralino with dominant SI 
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interaction and mass above 30 GeV have been discussed in ref. |34|, [35| , while the case 


37J. 


for an heavy neutrino of the fourth family has been considered in ref. 

II. WIMPs with mixed coupling in given model framework 

Since the 23 Na and 12 'I nuclei are sensitive to both SI and SD couplings - on the 
contrary e.g. of nai Ge and nat Si which are sensitive mainly to WIMPs with SI coupling 
(only 7.8 % is non-zero spin isotope in nat Ge and only 4.7% of 29 Si in ” at Si) - the analysis 
of the data has been extended considering the more general case of eq. (1) [[|l|]. This 
implies a WIMP having not only a spin-independent, but also a spin-dependent coupling 
different from zero, as it is also possible e.g. for the neutralino (see above). 

Following the usual procedure we have built the y log-likelihood function, which de¬ 
pends on the experimental data and on the theoretical expectations in the given model 
framework. Then, y has been minimized - properly accounting also for the physical con¬ 
straint set by the measured upper limit on recoils |22| - and parameters’ regions allowed 


at given confidence level have been obtained. In particular, the calculation has been per¬ 
formed by minimizing the y function with respect to the £asi, ^ctsd and mw parameters 
for each given 9 value. Here, ctsd is the point-like SD WIMP cross section on nucleon 
and tg6 is the ratio between the effective SD coupling constants on neutrons, a n , and on 
proton, a p ; therefore, 9 can assume values between 0 and 7r depending on the SD cou¬ 
pling. In the present framework the uncertainties on tq have been included; moreover, 
the uncertainties on the nuclear radius and the nuclear surface thickness parameter in 
the SI form factor, on the b parameter in the used SD form factor (see later) and on the 
measured quenching factors [Q of these detectors have also been considered 


Fig. 0 shows slices for some mw of the region allowed at 3 a C.L. in the (£ asi , £&sd-> 
mw) space for some 6 value. Only the case of four particular couplings are shown here 


for simplicity: i) 9 = 0 ( a n =0 and a p % 0 or 


>> a, 


J); ii) 9 = 7t/4 (a p = an); iii) 9 = 


7T/2 (a n % 0 and a p = 0 or 
The case a p = —a 


>> |a p |); iv) 9 = 2.435 rad (^ = -0.85, pure Z° coupling), 
is nearly similar to the case iv). 

Comments about other discussions of purely SD component as well as about the 
comparisons with other direct and indirect experiments can be found in ref. |H]. 

As already pointed out, when the SD contribution goes to zero (y axis in Fig. 0), an 
interval not compatible with zero is obtained for £a si■ Similarly, when the SI contribution 
goes to zero (x axis in Fig. 0), finite values for the SD cross section are obtained. Large 
regions are allowed for mixed configurations also for ^asi < 10” 5 pb and £,ctsd < 1 pb; 
only in the particular case of 9 = | (that is a p = 0 and a n % 0) £,(Jsd can increase up 
to ~ 10 pb, since the 23 Na and 12 'I nuclei have the proton as odd nucleon. Moreover, 
in ref. [[LL]] we have also pointed out that: i) finite values can be allowed for £,(Jsd even 
when £ctsi — 3 • 1CD 6 pb as in the region allowed in the pure SI scenario considered in the 
previous subsection; ii) regions not compatible with zero in the £,ctsd versus m w plane are 
allowed even when £crsi values much lower than those allowed in the dominant SI scenario 
previously summarized are considered; iii) minima of the y function with both £asi and 
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Figure 7: Slices for some mw of the region allowed at 3 a C.L. in the (£crsi, F.o'sd, mw) 
space for some 9 value in the model framework considered in ref. [[3TJ . Only four particular 
couplings are reported here for simplicity: i) 6 = 0; ii) 9 = 7t/4 iii) 6 = 7t/2; iv) 9 = 2.435 
rad. 


£ctsd different from zero are present for some mw and 9 pairs; the related confidence level 
ranges between ~ 3 a and ~ 4 cr [[H . 

In conclusion, this analysis has shown that the DAMA data of the four annual cycles, 
analysed in terms of WIMP annual modulation signature, can also be compatible with a 
mixed scenario where both £,<Jsi and £cr sd are different from zero. The pure SD and pure 
SI cases in the model framework considered here are implicitly given in Fig. [7] for the 
quoted mw and 9 values. 

Further investigations are in progress on these model dependent analyses to account for 
other known parameters uncertainties and for possible different model assumptions. As 
an example we recall that for the SD form factor an universal formulation is not possible 
since the internal degrees of the WIMP particle model (e.g. supersymmetry in case of 
neutralino) cannot be completely separated from the nuclear ones. In the calculations 
presented here we have adopted the SD form factors of ref. |37j estimated by considering 
the Nijmengen nucleon-nucleon potential. Other formulations are possible for SD form 
factors and can be considered with evident implications on the obtained allowed regions. 
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III. Inelastic Dark matter 


It has been recently suggested |38| the possibility that the observed annual modulation 
of the low energy rate could be induced by possible inelastic Dark Matter: relic particles 
that cannot scatter elastically off of nuclei. As discussed in ref. @], the inelastic Dark 


Matter could arise from a massive complex scalar split into two approximately degenerate 
real scalars or from a Dirac fermion split into two approximately degenerate Majorana 
fermions, namely x+ and X-, with a 5 mass splitting. In particular, a specific model 
featuring a real component of the sneutrino, in which the mass splitting naturally arises, 
has been given in ref. [[$|. It has been shown that for the x~ inelastic scattering on 
target nuclei a kinematical constraint exists which favours heavy nuclei (such as 12 'I) 
with respect to lighter ones (such as e.g. nat Ge) as target-detectors media. In fact, 
can only inelastically scatter by transitioning to x+ (slightly heavier state than x~) and 
this process can occur only if the X- velocity, v, is larger than: 


25 


Vthr ~ 


mwN 


( 2 ) 


where mwN is the WIMP-nucleus reduced mass and hereafter c — 1. This kinemati¬ 
cal constraint becomes increasingly severe as the nucleus mass, mjv, is decreased B- 
Moreover, this model scenario gives rise - with respect to the case of WIMP elastically 
scattering - to an enhanced modulated component, S m , with respect to the unmodulated 
one, So, and to largely different behaviours with energy for both S 0 and S m (both show 
a higher mean value) [[3^]. 

For the sake of completeness, we remind that - as stressed in ref. 
is suitable to reconcile in every case the DAMA result ||2B|. |7]. |2B 

B 


381 - this scenario 


CDMS-I claim 


however, it is worth to note that - as discussed e.g 


fTT| . pm with the 
in sect. Ill of 


ref. 


in reality the claim for contradiction made by CDMS-I was largely unjustified 


381 offers a further possible 


both for experimental and theoretical reasons. 

Anyhow, the proposed inelastic Dark Matter scenario 
model framework, which has also the merit to recover the sneutrino as a suitable WIMP 
candidate. 

Therefore, a dedicated energy and time correlation analysis of the DAMA experimental 
data has been carried out |32|] by considering the data collected during four annual cycles 
(statistics of 57986 kg x day ) f 


m 


301, |3T|, p2|. Here aspects other than the 


interaction type have been handled according to ref. |3T| , fixing in this way a given model 
framework. 

In this inelastic Dark Matter scenario an allowed volume in the space (^a p ;mw,5) is 
obtained |32]]. For simplicity, Fig. shows slices of such an allowed volume at some given 
WIMP masses (3 a C.L.). It can be noted that when mw the expected differential 

energy spectrum is trivially dependent on mw and in particular it is proportional to the 
ratio between £cr p and mw ; this particular case is summarized in the last plot of Fig. 
|l. The allowed regions reported there have been obtained by the superposition of those 
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Figure 8: Slices at fixed WIMP masses of the volume allowed at 3 a C.L. in the space 


See text. 


(£cr p , 5 , mw)'i the uncertainties on some of the used parameters have been included 32 


obtained when varying the values of the previously mentioned parameters according to 
their uncertainties. This also gives as a consequence that the cross section value at given 
5 can span there over several orders of magnitude. The upper border of each region is 
reached when Vthr approximates the maximum WIMP velocity in the Earth frame for 
each considered model framework (in particular, for each u 0 value). 

Note that each set of values (within those allowed by the associated uncertainties) for 
the previously mentioned parameters gives rise to a different expectation, thus to different 
best fit values. As an example we mention that when fixing the other parameters as in 
Ref. |3lH , the best-fit values for a WIMP mass of 70 GeV are: i) £ o v = 2.5 x 10~ 2 pb 
and S = 115 keV when v 0 = 170 km/s, ii) £cr p = 6.3 x 10 -4 pb and S = 122 keV when 


v 0 = 220 km/s. 

Finally, we note also here that significant enlargement of the given allowed regions 
should be expected when including complete effects of model (and related experimental 
and theoretical parameters) uncertainties. Moreover, possible different halo models can 
be also considered. 
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6.3 Conclusions 


The DAMA annual modulation data of four annual cycles |26|. [27], [28]. [Lj [31], [T2]| have 
been analysed by energy and time correlation analysis in terms of purely SI, purely SD, 
mixed SI/SD, “preferred” inelastic WIMP scattering model frameworks. 

To effectively discriminate among the different possible scenarios further investigations 
are in progress. In particular, the data of the 5 th and 6 th annual cycles are at hand, while 
the set-up is running to collect the data of a 7 th annual cycle. Moreover, the LIBRA 
(Large sodium Iodine Bulk for RAre processes) set-up is under construction to increase 
the experimental sensitivity. 


7 Towards LIBRA 

At present our main efforts are devoted to the realization of LIBRA (Large sodium Iodine 
Bulk for RAre processes in the DAMA experiment) consisting of ~ 250 kg of Nal(Tl). 
New radiopure detectors by chemical/physical purification of Nal and T1 powders as a 
result of a dedicated R&D with Crismatec have been realized. This will allow to increase 
the sensitivity of the experiment. 
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